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MODELING AND PARAMETER UNCERTAINTIES FOR 
AIRCRAFT FLIGHT CONTROL SYSTEM DESIGN 


by 

J. D. McDonnell, R. A. Berg, R. M. Heimbaugh, 
and C. A. Felton 

Douglas Aircraft Company 
McDonnell Douglas Corporation 
Long Beach, California 


SECTION 1 
INTRODUCTION 

Aircraft designs are continuing their evolution in structural and aerodynamic 
innovations that require artificial means of stability and dynamic augmenta- 
tion to obtain normal flying characteristics and acceptable structural 
margins. The motivation, of course, is to continue the reduction in direct 
operating costs for commercial aircraft and to continue to improve the mis- 
sion effectiveness/costs of ownership of military vehicles. The recent rapid 
rise in fuel costs has made commercial applications of some control- 
configured-vehicle (CCV) concepts desirable in the near term, perhaps on a 
derivative aircraft. Thus, it is not only imperative that dynamics of the 
aircraft in flight turn out as predicted, since stability margins will be 
lower (and in some cases, zero or negative), but that the flight control 
systems for those aircraft provide desired characteristics when the aircraft's 
actual dynamic description has some expected variability about its predicted 
behavior. 

The data presented in this report are intended to provide a start of a 
typical data base for use in flight controls work that is directed toward the 
kinds of systems noted above, i.e., those that are flight-critical or tending 
toward it. 

The data show typical historical variability in the steady aerodynamic and 
selected dynamic structural parameters of jet transports at one flight 
condition. Historical values of many dimensionless stability derivatives 


include original analytical estimates, wind tunnel updates, and final flight 
determined values. Structural dynamic parameters as originally estimated 
are shown, and the impact of configuration changes as a function of time are 
given. Final corrected values based on modal ground vibration test results 
are also given. 

The data base constructed for this report has been derived from many sources 
and is representative of many analytical and experimental techniques. The 
parameters themselves range from significant dynamic contributors to those 
with little influence on aircraft flying characteristics . The data, there- 
fore, indicate to the designer the historical uncertainties and variations 
that may be encountered in an advanced control system design and development 
program on a large subsonic transport aircraft. In a specific future design 
exercise, improvements in analytical or experimental techniques would have to 
be accounted for, as would greater uncertainties introduced as a result of 
new technology (and possibly less understood) devices and configurations. 

Section II of this report discusses the sources of aerodynamic, structural 
and weight data and defines the formats for data presentation. The methods 
used to derive the data are discussed as are the results themselves. The 
general data format objective has been to show the historical variations of 
the various parameters as a function of time. The historical variations of 
the dimensionless aerodynamic derivatives are given so that the effects of 
independent error sources on rigid body dynamics can be estimated regardless 
of the specific notation of the user (i.e., three degree-of-freedom linear 
small perturbation equations, six degree-of-freedom non-linear equations of 
motion, transfer functions, state space notation, etc.). 

Historical variations of structural dynamic characteristics are presented 
for the design evolution of large transport airplanes in terms of natural 
frequencies, the corresponding modal displacements of selected stations, and 
the corresponding generalized masses. The modal characteristics are dis- 
cussed with respect to the control designer as to value and as to limitations. 
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The important field of unsteady aerodynamics as related to structural res- 
ponse and flutter and control -configured vehicle design were not a part of 
this study but should not be neglected in future studies. 

Section III gives all of the actual data, while Section IV notes some possible 
applications of the data. Aerodynamic equations of motion and the definitions 
of the dimensional stability derivatives are given in the Appendix. 
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SECTION 2 

DISCUSSION OF DATA AND METHODOLOGY 


The data included in this report are from many sources and disciplines and 
cover time spans ranging up to four years. This section will discuss sources 
of data and presentation formats, and is divided into three parts where 
detailed discussions of aerodynamic parameter variations, and structural 
modal analysis are presented. Examples of data are given in this section to 
support the discussion, but the complete data set is given in Section III. 

AERODYNAMIC PARAMETERS 

The aerodynamic longitudinal and lateral-directional dimensionless deriva- 
tives given in Table 1 and 2 have been reconstructed as a function of time, 
where time in general goes back to initial configuration development for the 
particular airframe being considered. The relationship of these derivatives 
to aircraft dynamics is shown in the Appendix, where the aircraft equations 
of motion and their coefficients are defined in terms of the dimensionless 
derivatives. 

Detailed Examp 1 e of Stability Derivative Calculation 


An example of the type of detailed calculation used for the determination of 
a single stability derivative is presented here. The selected parameter is 
the fundamental derivative, (C[_ a )^ the airplane lift-curve slope, which can 
be calculated from the following expression: 


< C L 0 >A ' (Cl„) %* < c L„> ? 0 - { F H<CL„)^ /R (Cl,/ { { 1 -[(ft) * + ( 4 &)' 


Equation (1) 


where 


3C 


E ° L 

(Cl^a = f° r the entire airplane and includes static aeroelastic 


effects . 
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TABLE 1 

LONGITUDINAL AERODYNAMIC PARAMETER UNCERTAINTIES 


1 . 


2 . 


3 . 

4 . 

5 . 

6 . 

7 . 

8 . 
9 . 

TO. 

11 . 

12 . 

13 . 

14 . 

15 . 

16 . 
17 . 


Drag coefficient (Cp) 

Rate of change of drag coefficient with forward speed (C n ) 

. u u 

Rate of change of drag coefficient with angle of attack (C n ) 

u a 


Rate of change of drag coefficient with elevator deflection (Cn ) 

6 e 

Lift coefficient (C|_) 

Rate of change of lift coefficient with forward speed (Cl u ) 

Rate of change of lift coefficient with angle of attack (Ci ) 

Rate of change of lift coefficient with angle of attack rate (C. ) 

L a 

Rate of change of lift coefficient with pitch rate (C^) 

Rate of change of pitching moment coefficient with forward speed (C mu ) 

Rate of change of pitching moment coefficient with angle of attack 


( c m a ) 

Rate of change of pitching moment coefficient with angle of attack 

rate (C m .) 

a 

Rate of change of pitching moment coefficient with pitch rate (C m q) 
Rate of change of lift coefficient with elevator deflection (C^ ) 
Rate of change of pitching moment coefficient with elevator deflec- 
tion (C m & ) 
e 

Rate of change of drag coefficient with pitch rate (C^) 

Rate of change of drag coefficient with angle of attack rate (C D .) 
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1 . 
2 . 


5. 

6 . 

7. 

8 . 

9. 

10 . 

n. 

12 . 

13. 

14. 

15. 


17. 

18. 


TABLE 2 

LATERAL-DIRECTIONAL AERODYNAMIC PARAMETER UNCERTAINTIES 


Rate of 
Rate of 

(C n6 r ) 

Rate of 

\ ! 


change of yawing moment coefficient with sideslip angle (C n ) 
change of yawing moment coefficient with rudder deflection 

change of yawing moment coefficient with aileron deflection 


a 

4. Rate of change of yawing moment coefficient with spoiler deflection 

/ r- \ 


<V ' 

6SP • , 

Rate of change of yawing moment coefficient with yaw rate (C n ) 

Rate of change 


Rate of 
Rate of 

Rate of 

Rate of 

Rate of 

Rate of 

Rate of 

Rate of 

Rate of 


of yawing moment coefficient with roll rate (C np ) 

change of side force coefficient with sideslip angle (Cy g ) 
change 


change 

change 

change 

change 

change 

change 


of side force coefficient with rudder deflection (C ) 

^ 6 r 

of side force coefficient with aileron deflection (C ) 

y S, 

of side force coefficient with spoiler deflection (C ) 

y s 

of side force coefficient with yaw rate (Cy r ) sp 

of side force coefficient with roll rate (Cy p ) 
of rolling moment coefficient with sideslip angle (C^) 
of rolling moment coefficient with rudder deflection 


change of rolling moment coefficient with aileron deflection 


16. Rate of change of rolling moment coefficient with spoiler deflection 


(C W 

Rate of 
Rate of 


change of rolling moment coefficient with yaw rate (C^) 
change or rolling moment coefficient with roll rate (Cjj, p ) 
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(Cl o ) T0 


E/R _ 


*0 

Elastic 

\ 9a J 


( 3C l' 

Rigid 

_\ 9 a 



elastic effects to C. 

L a 

the tail . 


= the ratio of C, with aero- 

L a 

Tail -off 

of the rigid airplane, in the absence of 


(C L a ^TO 


^ r 

L for the rigid airplane in the absence of the tail. 
9a 


= factor that accounts for effect of aft fuselage bending on horizon- 
tal tail lift. 


(C L ) 


E/R = 
H 



Elastic 



Rigid 


elastic effects to C 
tal tail. 


= the ratio of C, with aero- 

L a 


Hori zontal 
Tai 1 

of rigid structure, of the isolated horizon- 


d 3 r 

(Cl ) = L horizontal tail lift curve slope of the rigid airplane 

“H 9 oij^ 


(^a) R = downwash gradient for the rigid airplane 

= increment to account for aeroelastic effects on downwash 
gradient 


The method of determination of the various components in the equation depends 
on the point in the design cycle at which the calculation is made. For this 
discussion, four stages in the design process will be considered: 
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1. Early preliminary design 

2. Later preliminary design, but before wind tunnel data are available 

3. After wind-tunnel data are available 

4. After flight-test data are available 

It should be noted that wind-tunnel tests are usually conducted prior to the 
official program start but for the flight condition selected in the present 
study these tests must be of the high-speed variety in order to correctly 
represent compressibility effects. 

In Stage 1 of the design effort the components of the equation which apply to 
the rigid ai-rplane are" found by relatively simple methods. For example, the 

D D 

derivatives (Cl a )y 0 and (Cl- ) can be determined through the use of the 
USAF DATCOM (Reference 1), from which Figure 1 is reproduced. The aero- 
elastic correction terms, (Cl^jo^ ( c L a ) E//R » and (A9e/9a) E , may be 
estimated based on values from previous similar production designs. 


During Stage 2 in the design more sophisticated methods are employed. Typical 
of these is the Weissinger lifting surface theory. Reference 2, which lends 
itself to the estimation of both rigid and elastic characteristics. 



FIGURE 1. SUBSONIC WING LIFT-CURVE SLOPE 
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Highspeed digital computers are a virtual necessity for the application of 
this and similar complex methods. Other aeroelastic calculations for wing 
and tail surfaces are made by means of the Hedman vortex lattice lifting 
surface theory for elastic wings, presented in Reference 3. 


The third stage of the design process benefits from the availability of wind- 
tunnel data. As previously noted, the tests must be appropriate for the 
flight condition of interest; i.e. high-speed wind tunnel data for cruise 
flight conditions for which compressibility effects are important, and high 
Reynolds number tests for high angle-of-attack flight conditions. Aside 
from the usual testing corrections that must be applied to wind tunnel data, 
other adjustments must be made to account for the fact that geometric dis- 
similarities often exist between the wind tunnel model and the configuration 
being analyzed. For the example derivative, (Cl )£, wind tunnel tests pro- 
vide the information for which the components (Cj_ )^q, (C, ) and (9e/3a) 

“ i H 


are derived. Ground tests to ascertain the airplane stiffness characteristics 
are conducted during this period to verify the estimated values. If any 
significant discrepancies exist, the aeroelastic corrections are reevaluated 
using the methods previously described. 


The final stage of the process follows the first flight and involves verifica- 
tion of the estimated values of the derivatives. Measurements made in flight 
test are obviously for the elastic aircraft, so the procedure is now reversed. 
In the case of the present example, the complete elastic derivative is 
measured, the elastic corrections are assumed accurate, and the elastic tail 


effectiveness F H (CL a )E /R (C|_i H ) 


can be measured, thus permitting the rigid 
airplane characteristics to be calculated. Any necessary adjustments are 
made to these rigid airplane derivatives. 


At this point some comments regarding stability derivatives and flying 
qualities may be appropriate. Prior to first flight in the development of 
an airplane, stability derivatives are carefully estimated as discussed in 
the preceding paragraphs. To the stability and control engineer the 
stability derivatives serve largely as the building blocks with which flying 
qualities are calculated. The flying qualities are then checked for com- 
10 



pliance against numerous criteria. Of course, the stability derivatives are 
also important to other engineering specialists, in particular, control sys- 
tem designers and others who utilize aircraft equations of motion in their 
work. 

Once the flight test program is underway, many of the airplane flying 
qualities can be determined directly without the necessity of knowing the 
values of specific derivatives. The situation is somewhat different than 
before flight test; it is generally the flying qualities information that is 
available and the estimated stability derivatives are checked using the 
test data. Flight test results, however, are obtained for discrete flight 
conditions and flying qualities analyses over the remainder of the flight 
envelope require knowledge of the stability derivatives. Therefore, the 
stability derivatives are established at the flight test points and inter- 
polations or fairings are made to develop continuous values. In this manner, 
a complete bank of aerodynamic data is compiled, from which the airplane 
flying qualities are generated for the entire flight envelope. 

Referring again to Equation (1) it is apparent that several of the terms are 
not unique to the expression for (Ci_ a l^; they are involved in the determina- 
tion of many of the other aerodynamic derivatives. Therefore, errors or 
uncertainties associated with these terms would cause a degree of correlation 
to exist between many of the derivatives. For example, the horizontal tail 
lift-curve slope CL a ^, the horizontal tail elastic-to-rigid correction factor 
(CL a ){!j » and the aft fuselage bending correction factor F H , are elements in 
the expressions for numerous longitudinal stability derivatives. In the data 
presented, the proliferation of an error in the Ci estimate is rather dra- 

a rl 

matic. The error (caused by an underestimate of Mach number effects on Cl^) 
is reflected in the estimates of seven separate stability derivatives: C[_ a , 
CL a’ C Lq» Cmq, CL fie , an< ^ ^ m 6e' P° tent i a ^ f° r this type °f error 

correlation exists with many of the other parameters, including the lateral - 
di recti onal variety. 

* 

So urces of Uncertainties 

There are a number of sources for error or uncertainty associated with the 
estimated values of the various aerodynamic parameters. In the very early 

11 


stages of design there is considerable likelihood that the airplane configu- 
ration will be altered before the first flight occurs. In some cases these 
modifications taks place well into the period between program commencement 
and first flight. Depending on the type of change made to the configuration, 
there can be a considerable impact on the estimated value of the aerodynamic 
parameters . 

Another uncertainty arises from inaccuracies that are inherent in analytical 
methods used to estimate the aerodynamic characteristics, including aero- 
elastic corrections. The degree of the uncertainty, when related as a 
percentage of the actual of final value, is magnified by the method of calcu- 
lating most derivatives as the sum of two or more increments. When two 
increments are of opposite sign and the magnitudes are similar, small errors 
in either increment results in a large error in the total value. For example, 
tail increments are often added to airplane tail -off increments to obtain 
total airplane derivatives, and depending on the particular derivative, the 
two increments may be of opposite sign. 

If the flight conditions of interest should change for any reason during the 
design process (e.g., modification of flight envelope), the coefficients of 
the equations of motion (dimensional stability derivatives) will change 
because, in general, they are functions of such flight parameters as Mach 
number, dynamic pressure, airspeed, altitude, and angle of attack. 

Then, of course, there are the ubiquitous computational errors that occasion- 
ally go undetected for significant periods of time. Obviously, it is not 
possible to predict the occurrence of uncertainties stemming from such causes. 

Discussion of Data 

The flight condition for which the subject data were compiled was deliberately 
chosen to be one which is particularly prone to aerodynamic modeling errors. 
This proclivity for error arises from the transonic flow conditions that 
exist at high-speed cruise and to the relatively high dynamic pressure which 
has a major effect on the aeroelastic corrections. Any small misjudgements 
in the effects of compressibility can be exaggerated because of the rapid 
variations with Mach number that occur with many aerodynamic characteristics 
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in the transonic speed range. A typical variation with Mach number of a 
primary stability derivative, C|_ a » the airplane lift-curve slope, is shown in 
Figure 2. The symbol denotes the approximate cruise flight condition of this 
s tudy . 

It should be observed that it is difficult to generalize on the aerodynamic 
parameter uncertainties. A derivative which is estimated with great accuracy 
on one airplane may not fare as well on the next. For example, uncertainties 
that are attributed to configuration changes certainly cannot be generalized. 
However, some generalizations associated with estimation accuracies can be 
offered. Table 3 a provides an indication of the relative impact and the 
estimation accuracy normally expected for each of the aerodynamic parameters. 

It should be emphasized that the categorizing of the parameters is approximate 
and will tend to vary with particular aircraft configurations and with flight 
conditions. It is emphasized that the estimation accuracy shown includes 
only the estimating methods normally used and does not consider uncertainties 
due to other causes. In Table 3 b the probable risk in flying qualities 
analyses and control system design is shown for each of the nine categories 
of Table 3 a. For example, a derivative with secondary impact that is esti- 
mated with only fair accuracy is no more or less likely to create problems 
than a primary derivative that normally is estimated with good accuracy. 


lIFT-CURVE 

SLOPE, 

(c lJto 



FIGURE 2. TYPICAL VARIATION OF LIFT-CURVE SLOPE WITH MACH NUMBER 
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TABLE 3 

PARAMETER ESTIMATION ACCURACY AND POTENTIAL 
ERROR IN FLYING QUALITIES ANALYSES 


IMPACT ON FLYING QUALITIES 

Primary 

Secondary 

Negligible 


(b) Potential Error In Flying 
Qualities Analyses: 

Large 

Moderate 
O Small 
Minimal 


SSI 











Historical values of each of the aerodynamic parameters are presented in 
Section 3. Most of the historical values are normalized by the final value 
of each parameter. However., in some cases in which the final level of the 
parameter is very small, the normalized value would be grossly distorted. 

In these cases the incremental difference between each value and the filial 
value is deemed more meaningful and is therefore presented. In the case of 
the derivative Cm^, the value depends on the selected moment reference cen- 
ter, which is another reason to present an incremental uncertainty. On each 
of the historical plots the abscissa is marked with the time of first flight 
(FF) and the FAA certification date (CERT), in addition to the authority to 
proceed (ATP, or program initiation) at the origin. Also noted on the plots 
is the time of the first complete dynamics analysis (D) of the airplane. 

This date is significant because it represents the first firm requirement 
for some of the parameters. Included with each graphical presentation are 
brief comments regarding the usual impact of the particular parameter on the 
augmented airplane dynamics. Also included are comments relating to the 
final value uncertainty. A typical presentation is repeated here for conven- 
ience as Figure 3. 

In viewing the data, certain of the parameters reveal a relatively large 
degree of uncertainty early in the design program. In most of the cases, 
this is attributable to configuration changes being made as the design is 



PRIMARY AERODYNAMIC DERIVATIVE 
CONTRIBUTES TO DAMPING OF SHORT- 
PERIOD MODE 

ESTABLISHES LEVEL OF ACCELERATION 
SENSITIVITY, n/o (NUMERATOR 
PROPERTY OF olb e TRANSFER FUNCTION) 
FINAL VALUE VERY CAREFULLY CHECKED 
IN FLIGHT TEST BECAUSE OF PERFORM- 
ANCE GUARANTEES 


FIGURE 3. TYPICAL DATA FORMAT; C. „/C, „ - 

La La FINAL 
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finalized. High levels of uncertainty near the outset of the program in 
some cases can be traced to the lack of early aeroelastic corrections; that 
is, the available estimates of the parameter are for the rigid airplane. Had 
there been a firm requirement for these parameters at this time some approx- 
imate elasticity corrections would have been applied. Similarly, compres- 
sibility correction were not applied during the initial design stages to some 
of the less significant parameters. 

It will be noted that, in general, the data points do not tend to coincide 
chronologically for the various parameters. The point in time at which a 
derivative is estimated is affected by several factors, among which are the 
importance of the derivative and the difficulty of estimation. When a con- 
figuration change occurs during the design the more important derivatives are 
recalculated first. In fact, configuration changes are not made until the 
impact on certain parameters and aircraft flying qualities is ascertained. 

Of the aerodynamic parameters presented in Section 3 all except two are 
partial derivatives. These two are the airplane lift and drag coefficients 
which are given because they appear in the coefficients of the equations of 
motion, sometimes referred to as the dimensional stability derivatives (see 
Appendix). The airplane lift coefficient is an independent variable in 
specifying a particular flight condition so no uncertainty can be assigned. 

As discussed earlier, however, if the flight condition of interest should 
change for any reason, not only the lift coefficient, but the entire set of 
equations of motion will be altered. 

The aircraft drag coefficient is somewhat unique in this discussion because 
it is the basic aerodynamic design parameter for cruise flight. The initial 
value shown in Section 3 represents a guarantee and normally is not revised 
until flight test data are available. Every effort is made to achieve this 
design goal, including configuration changes which can range from minor 
refinements in the design to major drag reduction programs. 

The uncertainties associated withthe final values of the various aerodynamic 
parameters are difficult to assess when traditional methods are used in the 
analysis of flight test data, as in the subject case. The measurement 
accuracy of certain of the parameters receives great attention because of the 
relationship to performance guarantees. An example is presented in Figure 4 
16 




where drag coefficient, Cd» has been plotted versus Mach number at various 
levels of lift coefficient. Cl. The scatter in the data can be regarded as 
an indicator of the uncertainty in the final value; for the subject case, the 
standard deviation of Cg is 2 percent. The total error arises from several 
factors, the major one being the engine thrust measurement error. 

There are other parameters which are also carefully checked against flight 
test data but are not normally analyzed with respect to measurement error. 

An example is presented in Figure 5 where the horizontal stabilizer effec- 
tiveness, Cnrj H , is shown as a function of Mach number. This parameter is 
simply related to the elevator control power term, C m(S , so that the measure- 
ment accuracy of the two derivatives is comparable. In Figure 5, only data 
points near the subject flight condition are shown, and the calculated stan- 
dard deviation over this Mach number range is 4 percent. An insight into the 
potential measurement error of C m -j^ can be gained by considering the expres- 
sion used to calculate the individual points: 
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FIGURE 5. TYPICAL FLIGHT TEST DATA; STABILIZER EFFECTIVENESS 


where the incremental values represent the differences between two trimmed 
centers of gravity. Errors can arise from several sources: the measurement 
of the stabilizer angle; the measurement of the eg location; the measure- 
ment of airplane weight and dynamic pressure that determine the lift 
coefficient; and the ability of the pilot to precisely trim the aircraft. 

The level of measurement error for H is considered typical of a number of 
static derivatives. 

In the case of the dynamic derivatives, the verification is normally accom- 
plished by a comparison of response characteristics; i.e. frequencies, 
damping ratios, time to half or double amplitude, etc. If the match between 
the estimated and the measured characteristics is reasonable, the estimated 
derivatives are presumed to be correct; if not, adjustments are made to the 
derivatives in order to achieve a satisfactory match of the aircraft dynamics. 
For the most part, it is extremely difficult, if not impossible, to isolate 
the undertainties or errors associated with the measurement of individual 
dynamic derivatives. 
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WEIGHT AND INERTIAL PARAMETERS 


Aircraft weight and inertias are key quantities i r. the aircraft equations of 
motion and have a history of variability during the aircraft's configuration 
development. Although the aircraft can be weighed quite accurately prior 
to first flight (to within +_ 0.2 percent), the manufacturer's empty weight 
(MEW) and gross weight can change significantly after program initiation. 
Manufacturer's empty weight changes occur both because of changes in the 
specification (increase in range, say) and changes in material (engine weight 
change). If the designer's concern is control system dynamics, he will be 
more interested in gross weights than MEW's. Both are given, however, so 
that the user of these data will have some appreciation for the kinds of fac- 
tors influencing his control system design. 

The growth of MEW over the four-year period is shown in Figures 6 and 7, 
where it will be noted that the initial value was 7.7 percent less than the 
final value. The major jumps in MEW are accounted for in Table 4. Of the 
7.7 percent difference, specification growth accounted for 4.7 percent and 
non-specification growth 2.0 percent. Figure 8 gives the changes in other 
design weights for the same period of time. Note from Figure 6 that the 
original maximum takeoff gross weight (MTOGW) was 90 percent of the final 
val ue. 

Coincident with the weight increases are changes in moments of inertia. 

These are summarized in Table 5. It will be noted that the increase in gross 
weights can be utilized several ways. The operation can keep the fuel load 
constant and increase payload, in which case the pitching moment of inertia 
Iy, shows a rather large increase (Table 5, Item 2). When the payload is held 
constant and fuel load is increased, the roll moment of inertia, 1^, shows 
a large increase, as shown in Table 5, Item 3. 

STRUCTURAL PARAMETERS 

Modal Vibration Analysis 

The approach for determining the uncertainties associated with structural 
dynamic considerations arising from structural flexibility was to track his- 
torically the variations in the normal modes of vibration as described by 
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SEE TABLE 4 FOR DESIGN WEIGHTS AND BRIEF SUMMARY OF INDICATED CHANGES 










TABLE 4 

MAJOR WEIGHT CHANGE SUMMARY 

(see Figure 6) 


A. Revised estimate for wing bending material, 
increased allowance for insulation and 
interior panels plus miscellaneous changes. 

B. Increased engine weight, revised estimate 
for wing and landing gear to provide for an 
increase in takeoff gross weight plus 
miscellaneous changes. 

C. Interior design changes, 0.102m (4-inch) 
fuselage stretch, addition of aft engine 
maintenance platform plus revised estimates 
for fuselage, wing and tail. 

D. Customer requested interior changes plus 
revised estimates for miscellaneous struc- 
tural and subsystem items. 

E. Increased engine weight, revised estimate 
for wing and landing gear to provide for an 
increase in takeoff gross weight plus mis- 
cellaneous changes. 

F. Incorporated more representative weights 
for galleys and passenger seats. 


frequencies, modal displacements, and the resulting generalized masses. The 
method for computing the modes was begun by the usual lumping of the fuselage, 
empennage, and wing into a finite number of mass and spring elements. As the 
aspect ratio of the selected airplane was large, it was justifiable to 
represent the structure as interconnected slender beams. The root or 
connecting springs were computed from a finite element representation using 
a redundant force procedure. 

The time period considered in this study covered 35 months starting with the 
formal authority to proceed (ATP) and ending with the ground vibration test 
(GVT) 1 month prior to first flight. Preceding this 35-month period there 
was approximately 6 months of preliminary design activity. This preliminary 
design phase was not covered by this study because structural uncertainties 
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TABLE 5 


MOMENT OF INERTIA HISTORY 

Inertia Value at Program Start 
Assumed Growth Final Inertia Value 

1. Manufacturers Empty Weight (MEW) 


Growth of 8.8 percent 


V ! X Final 

0.95 

V ! Y Final 

0.92 

V*Z Final 

0.93 

Maximum Takeoff Gross Weight 
(MTOGW) Growth of 11 percent 
and Constant Fuel: 


V T X Final 

0.96 

V*Y Final 

0.88 

l l /l l Final 

0.92 

Maximum Takeoff Gross Weight 
(MTOGW) Growth of 11 percent 
and Constant Payload: 


V ! X Final 

0.85 

Final 

0.93 

l l ll l Final 

0.89 


Final values are those calculated at first flight (31 months 
after program initiation). The product of inertia 
(I xz ) does not change significantly. 
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would have been distorted by the large configuration changes resulting from 
changes in design specifications that were made during this phase. 

During the evolution of the design from ATP to GVT eleven changes were made 
which were judged of sufficient significance to warrant computation of new 
modes. During that time, the modes were used primarily for flutter analyses 
with application to control system effects being made less frequently. For 
this study the modes were recomputed and the results compared to the final or 
certification results. The final results were computed based upon corrections 
as made from GVT results. 

In particular, the dynamic representation of the airplane consisted of a 
masswise representation composed of 50 lumped bays. Each bay was capable of 
six "rigid body" degrees of freedom, three translations and three rotations. 
The entire mass of the airplane was lumped into these bays by calculating the 
mass properties of each bay about a selected reference station in the bay. 

Mass and stiffness symmetry was assumed allowing an analysis of half the 
airplane. Figure 9 shows the mass bay distribution used in the modal 
calculations . 

Theoretical cantilevered vibration modes were calculated using the stiffness 
and mass data of the respective bays for the major components - wing, hori- 
zontal stabilizer, fuselage and vertical stabilizer. All modes were 
calculated using structural influence coefficients (SICs) which relate con- 
trol point static deflections to applied forces. Complex joint structures 
such as wing-fuselage, fin-fuselage, and engine pylons, used structural 
influence coefficients obtained from finite element redundant force type 
analyses. Figure 10 illustrates the stiffness representation used in the 
vibration analyses. 

Rigid body modes were used in conjunction with the above component modes to 
release the airplane and generate free-free orthogonal modes into symmetric 
and antisymmetric sets. In all cases, sufficient component modes were used 
in the free-free modal calculations to assure modal convergence of the air- 
plane modes of concern, from zero to 10 hertz. All modal calculations were 
performed using a digital computer program. 
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It should be noted that for this study all the vibration analyses were per- 
formed using the final structural representation. That is, the early modes 
were reconstructed using the respective stiffness and mass data, but with the 
final 50 bay structural representation. In actual practice, the earlier 
vibration analyses were performed using a much coarser structural . represen- 
tation. Therefore, the uncertainties shown for the early vibration analyses 
for frequencies greater than 3 or 4 hertz may be greater than these studies 
indicate. 

As noted, the structural dynamic uncertainties were limited to those result- 
ing from mass and stiffness variations only. In the practical consideration 
of structural dynamic effects, whether for the usual application to flutter 
or gust loads or for control -continued vehicle application, the predictability 
of the unsteady aerodynamic terms is as important as the predictability of 
the mass and stiffness dependent modes of vibration. Especially important 
are the unsteady aerodynamics of control surfaces which must be included in 
the analysis of all aircraft with highly responsive active control systems in 
order to assure that the required flutter speed and gust load margins are 
met and that flutter and excessive gust and oscillatory loads are avoided 
under any likely control system failure or malfunction. 

Structural Dynamic Parameters Investigated 

During the design phase, the airplane design data are continually changing 
due to configuration changes, revisions based upon test data, or refinements 
to existing data based upon detailed analyses. Many of these changes are 
minor and do not significantly affect the airplane vibration characteristics. 
Therefore, only those parameters which significantly impacted the airplane 
vibration modes were investigated. These parameters were as follows: 

• Wing stiffness 

• Fuselage stiffness 

• Wing engine pylon stiffness 

• Aft engine pylon stiffness 

• Wing engine weight 

• Wing engine eg location 

• Aft engine weight 
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The chronology and relative magnitude of these parameter changes are shown 
in Table 6. Symmetric and antisymmetric vibration modes were calculated for 
both empty and full fuel conditions. These two fuel conditions result in the 
highest and lowest modal frequencies of the system and therefore are typical 
of the modal uncertainties which can be expected over the entire fuel range. 
For both fuel conditions the payload was held constant. These payload and 
fuel conditions match configurations used during the airplane ground vibra- 
tion test. Thus, a direct comparison of theoretical and experimental modal 
vibration data was easily made which facilitated corrections to the 
theoretical data. 

TABLE 6 

DESIGN PARAMETER SUMMARY 


DESIGN PARAMETER CHANGES 


TIME 

MON 


0 

2 

4 

5 

6 
7 

10 

12 

14 

16 

20 

24 

35 


WING 

STIFFNESS 


BASIC 

+30% BEND 
+ 13% TOR 


+5% BEND 
+ 13% TORS 


-2% BEND 
+6% TORS 

+1% BEND 
+1% TORS 

-2% BEND 
-2% TORS 

+5% TORS 


FUSELAGE 

STIFFNESS 


BASIC 


-2% TORS 


-2% BEND 
-2% TORS 


WING 

PYLON 

STIFFNESS 


BASIC 


+64% 


+37% 


+ 6 % 


AFT 

PYLON 

STIFFNESS 


BASIC 


T 

-35% 


t 

• 32 % 


WING ENG 
WEIGHT 
AND CG 


BASIC 


0.81m (32 IN.) AFT 
CG SHIFT 


24% WEIGHT 
INCREASE 


AFT 

ENGINE 

WEIGHT 


BASIC 


6% WEIGHT 
INCREASE 


NOTE: THE ABOVE CHANGES IN EACH PARAMETER ARE CUMULATIVE 
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Structural Dynamic Data 


The results of the modal vibration uncertainty study are presented as plots 
of modal frequency, generalized mass and displacement versus time. The 
historical evolution of these parameters is shown for a selected subset of 
the airplane vibration modes which are important in both control system 
design, gust loads due to turbulence, and flutter calculation, all with or 
without control -configured vehicle technology. The chronology of changes to 
the input data is shown on each plot for easy reference. The mode names, 
such as first wing bending, are chosen so as to describe the predominant 
motion of each mode although, of course, the motion involves the entire air- 
plane structure. Each parameter is given as a ratio of the value at a given - 
time to the final value. The final value of each parameter is taken from a 
theoretical set of modes which are in agreement with ground vibration test 
resul ts . 

The generalized mass and deflection plots for each mode were calculated using 
a consistent normalization point for each particular mode. For instance, 
first wing bending was normalized to unit vertical deflection at the wing tip 
bay. The modal displacements are presented at selected stations which repre- 
sent a variety of possible sensor locations which might be selected to 
implement various control systems. The following locations and deflections 
were selected: 

Symmetric Analysis 

1. Fuselage nose pitch angle (<*) 

2. Fuselage eg vertical deflection (h) 

3. Wing tip vertical deflection (h) 

4. Wing tip pitch angle (<*) 

5. Horizontal stabilizer tip vertical deflection (h) 

6. Horizontal stabilizer tip pitch angle [a) 

Antisymmetric Analysis 

1. Fuselage nose roll angle (e) 

2. Fuselage nose yaw angle ( ^ ) 

3. Wing tip vertical deflection (h) 
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4. Wing tip pitch angle (a) 

5. Horizontal stabilizer tip vertical deflection (h) 

6. Horizontal stabilizer tip pitch angle ( a ) 

7. Vertical stabilizer tip lateral deflection («,) 

8. Vertical stabilizer tip yaw angle ( <p ) 

The modal data are presented in Section 3. The type of data shown for each 
airplane configuration is summarized in Section 3. 

Discussion of Results 

General . An examination of the modal frequency, generalized mass and 
displacement plots show, in general, that these uncertainties can be divided 
into two distinct phases. The first phase occurs during the early design 
stages, approximately 0 to 10 months, where the modal parameters exhibit 
large variations. These variations reflect major configuration changes such 
as engine weight, engine location, maximum gross weight and the introduction 
of stiffness constraints to satisfy dynamic conditions such as landing, gust 
and flutter requirements. In addition, data generated during this phase are 
based upon a coarse idealization of the structure. 

The second phase occurs from approximately 10 to 24 months, during which time 
the modal parameters, in general, exhibit only small changes. This reflects 
the fact that the airplane configuration has been determined and all major 
design constraints have been introduced. Changes to the design data are 
minor modifications based upon additional detailed structural analyses, which 
produce small changes in the airplane modal characteristics. 

But however detailed the analysis, uncertainties still exist. This can be 
seen by comparing the modal parameter values at the 24th month with the 
certification values. The 24th-month points represent the final predicted 
theoretical modes before verification by ground vibration test. The theore- 
tical modes are a result of a highly detailed structural analysis of a fixed 
configuration. Therefore, the difference between the 24th month and certi- 
fication values represent uncertainties in analysis techniques in idealizing 
the airplane stiffness and mass properties. The idealizations must reduce 
the continuous airplane structure to a finite number of elements, making 
simplification inevitable. Also* the. number of elements used to represent 
the airplane is limited by practical computational and economic considerations 



An examination of the modal frequency, generalized mass and displacement 
plots does not show any discernable trend due to fuel condition or symmetry. 
The correlation is equally good for zero and full fuel and for symmetric and 
antisymmetric conditions. 

Modal Frequencies . The modal frequency plots. Section 3, show that extreme 
variations of -30 to +84 percent may occur for some modes during the first 
phase. However, the majority of the modal frequencies appear to be in the 
+20 percent range. During the second phase, the modal frequencies show 
variations of only 3 to 4 percent on the average. The final theoretical 
modes (24 month point) differ from the certification modes by 4 to 5 percent 
on the average, although maximum differences reach +13 to -20 percent for a 
small percentage of the modes. An average of 5 percent difference in modal 
frequency is extremely good correlation considering the complexity of the 
structure being analyzed. 

Generalized Mass and Displacement . The modal generalized mass and displace- 
ment plots. Section 3, show much the same trend as do the frequency plots. 
However, the magnitudes of the variations are much more extreme. It appears, 
in general, that changes in design data produce larger variations in 
generalized mass and deflections than in modal frequency. 

However, care must be exercised in interpreting the generalized mass and 
deflection data. It is true that the data show the uncertainty of each para- 
meter as a ratio of its value at a given time to its final value. But the 
actual value of the parameter may be small. Therefore, seemingly large 
variations in its value may still be insignificant as they affect the air- 
plane response characteristics and hence the control system design. In 
addition, the ratio of any particular parameter, such as fuselage nose pitch 
angle, varies for each mode of a given configuration. Note that the para- 
meter ratio is greater than unity in some modes and less than unity in others 
in which there is usually a bias in one direction or the other. 

Evaluation of Results . The importance of the structural dynamic behavior of 
the aircraft to the control system design depends upon the intended function 
of the control system and its gain/phase properties. For low frequency con- 
trol systems such as conventional autopilot and yaw damper systems, the 
effect of the aircraft stiffness distribution on the steady aerodynamic 
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derivatives is of significance. For higher frequency control systems, 
particularly those required for gust load and flutter suppression, the 
stiffness and mass distributions become more significant and hence the 
vibration modes of the aircraft become meaningful measures of potential 
control system interaction. Of course, in any control loop with significant 
gain and phase shift in the elastic modal frequency range knowledge of the 
vibration characteristics of the vehicle is desirable. Use of mode shapes 
to determine nodal and anti-nodal points for sensor locations is desirable 
for passive gain stabilization procedures, but for large flexible aircraft 
significant shifts in node lines can occur with fuel and payload distribution 
changes. 

The use of natural frequencies to assess the change in the vibration modes 
during the design cycle is a useful (but not unique) measure of potential 
effects on the control system design. The natural frequencies may be iden- 
tified by the aircraft components which are the most directly involved in the 
mode and some subjective evaluation will be made by the engineer as to the 
significance to any given control system. For example, a low frequency 
fuselage bending mode may have a very direct bearing on the control system 
design for a high gain system and a low frequency is an indication that 
significant data error from forward mounted sensors may occur in turbulence 
conditions. 

Knowledge of amplitudes at various stations on the aircraft and the general- 
ized masses is considerably less useful. The control system designer faced 
with the task of designing a system for response or stability in the elastic 
frequency range or a system with significant gain and phase variation in this 
frequency range must use these data in his analyses along with the appropriate 
aerodynamic functions. 

A review of the modal parameters variations through the design cycle presented 
herein indicates a reasonable convergence with time on the final design 
values for the modal frequencies. The generalized mass and modal deflections 
for various modes and locations on the aircraft, however, indicate large 
variations and an apparent lack of continuity. This latter is not surprising 
in view of the fact that during the design cycle data updates resulting 
from structural modifications and improved data (resulting from both quality 
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and idealization modifications) are generally made for one or two aircraft 
components at a time. For example, a shift in wing engine eg with respect 
to the wing may not significantly effect the basic fuselage and empennage 
modes but can cause large variations in generalized mass of modes involving 
significant engine motion. The frequency variation however will show signi- 
ficantly less variation. 

The importance to the control system design from the above change will be 
insignificant in the basic low frequency autopilot design but of major signi- 
ficance for a flutter suppression system. 

A more significant measure of the structural dynamic parameter uncertainties 
may be found by observing the variation of the aircraft open loop transfer 
function for specific control surface inputs at specific potential sensor 
locations. For such response analyses to be valid in the elastic frequency 
ranges, they must use the unsteady aerodynamic functions which more realis- 
tically represent the complex aerodynamic forces which occur in these 
frequency ranges. 

To illustrate the above, a frequency response analyses was performed to 
determine if the modal data herein would correlate with the aircraft transfer 
functions for two points in the design cycle. The transfer functions were 
evaluated for 0 to 10 Hertz for a plus-or-minus 1-degree elevator oscillation 
using unsteady aerodynamics for Mach 0.88 and flight at 7,315 meters (24,000 
feet). Table 7 shows a summary of results for two dynamic aeroelastic modes, 
the aircraft short-period mode and the Wing Engine Pitch/Fuselage Bending 
mode. The input modal displacement and generalized mass ratios and output 
acceleration ratios for these modes are shown for several aircraft locations. 
Free-free aircraft modes were used in the analysis. 

It is apparent from these summary data that while the generalized mass and 
normalized modal displacements show large variations, the acceleration 
responses for the selected locations are similar and indicate relatively 
small variations. These response data are of significant use to the control 
system designer using the elevator for force generation and any of the 
selected location response parameters for sensor locations. 
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TABLE 7 

RESPONSE ANALYSIS SUMMARY 


SYMMETRIC ZERO FUEL 

MACH 0.88 731 SM (24,000 FT) 


RESPONSE 

PARAMETER 

3.23 Hz 

MODE INPUT MODEL 
DISPLACEMENT 
RATIO 

RESPONSE ACCEL. RATIO 
< *24 / *e ) 

0.45 Hz 

SHORT PERIOD MODE 

3.23 Hz 

FUSELAGE BENDING MODE 

FUSELAGE NOSE a 

0.80 

1.02 

1.18 

FUSELAGE CG h 

3.21 

1.02 

1.12 

WING TIP h 

0.94 

1.02 

1.04 

WING TIP a 

0.44 

1.03 

0.85 

HORIZONTAL 




STABILIZER h 

1.00 

1.02 

1.27 


• ABOVE RATIOS ARE FOR 24TH MONTH DATA VERSUS CERTIFICATION DATA 

• (GM 24 /GM c ) = 0.51 FOR THE 3.23 HERTZ FUSELAGE BENDING MODE 


For any given aircraft structure the structural dynamic characteristics vary 
significantly over the normal fuel and payload range. As a matter of interest 
the symmetrical modal frequency variation from full fuel to zero fuel was 
compared to the maximum frequency variation for zero fuel over the design 
period from ATP to GVT. A similar comparison was made for the full fuel 
frequencies over the design period. These comparisons are shown in Table 8 
and indicate that the modal frequencies vary with fuel loading by the same 
level of magnitude as the uncertainty variations and in some cases more. 
Payload variations, which have not been considered, would show even a larger 
range in the frequencies of the final design. 
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TABLE 8 


COMPARISON OF FUEL EFFECTS TO UNCERTAINTIES EFFECT 
SYMMETRIC MODAL FREQUENCIES 


MODE DESCRIPTION 

f FULL FUEL 
f ZERO FUEL 

UNCERTAINTIES RATIO 
f MAX /f MIN 

ZERO FUEL 

FULL FUEL 

FIRST WING BENDING 

0.63 

1.26 

1.16 

WING ENGINE YAW 

0.98 

1.20 

1.16 

FUSELAGE BENDING 

0.90 

1.13 

1.09 

HORIZONTAL STABILIZER BENDING/AFT ENG PITCH 

1.01 

1.72 

1.08 

WING INNER PANEL TORSION 

0.90 

1.39 

1.35 

WING ENGINE ROLL 

1.01 

1.16 

1.15 

SECOND WING BENDING 

0.71 

1.30 

1.09 

WING FORE AND AFT BENDING 

0.56 

1.27 

1.20 

WING TORSION/ENGINE PITCH 

0.94 

1.13 

1.17 




SECTION 3 
DATA 


Aerodynamic, Inertial, and structural data are presented in this section. 

The data are discussed in detail in Section 2 so that this section provides 
as concise a compendium as possible. The parameters are defined in the list 
of symbols at the beginning of this report. The equations of motion and 
dimensional stability derivative definitions are given in the Appendix. The 
data are organized as follows: 


Longitudinal Aerodynamic Parameters : 


Figure 11 Historical Uncertainties of C, , 

L a 

Figure 12 Historical Uncertainties of C. , 

m Of 

Figure 13 Historical Uncertainties of C Q , 

u u 

Figure 14 Historical Uncertainties of C 

m S e 




e 


Lateral -Directional Aerodynamic Parameters: 


Figure 15 Historical Uncertainties of 
Figure 16 Historical Uncertainties of 
Figure 17 Historical Uncertainties of 

Figure 18 Historical Uncertainties of 

Figure 19 Historical Uncertainties of 

Inertial Parameters : 

Figure 20 Manufacturer's Empty Weight 
Figure 21 MEW Growth History - Cabin, 
Table 9 Major Weight Change Summary 
Table 10 Moment of Inertia History 
Figure 22 Design Weights Histories 



Growth History 

Wing, Empennage, Forward Fuselage 



Structural Dynamic Parameters : 

Table 11 Symnetrlc Modal Analysis Summary 
Table 12 Antisymmetric Modal Analysis Summary 
Figures 23 through 26 Mode Shapes and Mode Lines 

Figures 27 through 60 Evolution of Frequency Ratio (see Tables 11 and 12 
for index) 

Figures 61 through 90 Evolution of Generalized Mass Ratio (see Tables 11 and 
12 for index) 

Figures 91 through 196 Evolution of Displacement Ratio (see Tables 11 and 12 
for index) 
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FIGURE 11. HISTORICAL UNCERTAINTIES OF LONGITUDINAL AERODYNAMIC PARAMETERS 
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PRIMARY DIRECTIONAL CONTROL 
DERIVATIVE 

FINAL VALUE CHECKED BY FLIGHT 
TEST RESULTS; NORMALLY NO 
ASSESSMENT OF ACCURACY MADE 


FIGURE 17. HISTORICAL UNCERTAINTIES OF LATERAL-DIRECTIONAL AERODYNAMIC PARAMETERS 
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SECONDARY CONTROL DERIVATIVE 
FINAL VALUE NOT CHECKED IN 
FLIGHT TEST; IMPACT TOO SMALL 
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NEGLIGIBLE, NOT PLOTTED 


NEGLIGIBLE IN MANY CASES 
FINAL VALUE NOT CHECKED IN 
FLIGHT TEST; IMPACT TOO SMALL 
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SECONDARY CONTROL DERIVATIVE 
AFFECTSTURN COORDINATION 
CHARACTERISTICS 
FINAL VALUE NOT CHECKED IN 
FLIGHT TEST; IMPACT TOO SMALL 
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PRIMARY LATERAL CONTROL 
DERIVATIVE 

FINAL VALUE CHECKED BY FLIGHT 
TEST RESULTS; NORMALLY NO 
ASSESSMENT OF ACCURACY MADE 


NEGLIGIBLE, NOT PLOTTED 


sp 

NEGLIGIBLE IN MANY CASES 
FINAL VALUE NOT CHECKED IN 
FLIGHT TEST; IMPACT TOO SMALL 


FIGURE 18. HISTORICAL UNCERTAINTIES OF LATERAL-DIRECTIONAL AERODYNAMIC PARAMETERS 
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SECONDARY CONTROL DERIVATIVE 
AFFECTS TURN COORDINATION 
CHARACTERISTICS 
FINAL VALUE NOT CHECKED IN 
FLIGHT TEST; IMPACT TOO SMALL 
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DERIVATIVE 

FINAL VALUE CHECKED BY FLIGHT 
TEST RESULTS; NORMALLY NO 
ASSESSMENT OF ACCURACY MADE 


FIGURE 19. HISTORICAL UNCERTAINTIES OF LATERAL-DIRECTIONAL AERODYNAMIC PARAMETERS 
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TABLE 9 

MAJOR WEIGHT CHANGE SUMMARY 

(see Figure 6) 

A. Revised estimate for wing bending material, 
increased allowance for insulation and 
interior panels plus miscellaneous changes. 

B. Increased engine weight, revised estimate 
for wing and landing gear to provide for an 
increase in takeoff gross weight plus 
miscellaneous changes. 

C. Interior design changes, 0.102m (4-inch) 
fuselage stretch, addition of aft engine 
maintenance platform plus revised estimates 
for fuselage, wing and tail. 

D. Customer requested interior changes plus 
revised estimates for miscellaneous struc- 
tural and subsystem items. 

E. Increased engine weight, revised estimate 
for wing and landing gear to provide for an 
increase in takeoff gross weight plus mis- 
cellaneous changes. 

F. Incorporated more representative weights 
for galleys and passenger seats. 
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TABLE 10 

MOMENT OF INERTIA HISTORY 


Assumed Growth 

1. Manufacturers Empty Weight (MEW) 
Growth of 8.8 percent 

V ! X Final 

V ! Y Final 

l Z /l l Final 

2. Maximum Takeoff Gross Weight 
(MTOGW) Growth of 11 percent 
and Constant Fuel: 

V ! X Final 

V*Y Final 

h ,l l Final 

3. Maximum Takeoff Gross Weight 
(MTOGW) Growth of 11 percent 
and Constant Payload: 

V ! X Final 

V ! Y Final 

H'h Final 


Inertia Value at Program Start 
^ Final Inertia Value 

0.95 

0.92 

0.93 

0.96 

0.88 

0.92 

0.85 

0.93 

0.89 


Final values are those calculated at first flight (31 months 
after program initiation). The product of inertia 
(1^) does not change significantly. 
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FIGURE 22. DESIGN WEIGHT HISTORIES 




TABLE 11 


SYMMETRIC MODAL ANALYSES SUMMARY 


FIGURE NUMBER OF TIME HISTORY 



ZERO FUEL 

FULL FUEL 

MODE DESCRIPTION 

FREQ 

GEN 

MASS 

DEFL 

FREQ 

GEN 

MASS 

DEFL 

FIRST WING BENDING 

27 

" 

81 

91-96 

36 

69 • 

115-120 

WING ENGINE YAW 

28 

62 

- 

37 

70 

- 

WING ENGINE PITCH/ 
FUSELAGE BENDING 

29 

63 

97-102 

38 

71 

121-126 

HORIZONTALSTABILIZER BENDING/ 
AFT ENGINE PITCH 

30 

64 

- 

39 

72 

- 

WING INNER PANEL TORSION 

31 

65 

103-108 

40 

73 

127 131 

WING ENGINE ROLL 

32 

- 

- 

43 

- 

- 

SECOND WING BENDING 

33 

66 

109-114 

42 

75 

132-135 

WING FORE AND AFT BENDING 

34 

67 

- 

41 

74 

- 

WING TORSION/ENGINE PITCH 

35 

68 

- 

45 

- 

- 

THIRD WING BENDING/ 
AFT ENGINE PITCH 

- 

- 

- 

44 

'76 

- 


TABLE 12 

ANTISYMMETRIC MODAL ANALYSES SUMMARY 




FIGURE NUMBER OF TIME HISTORY 




ZERO FUEL 


FULL FUEL 

MODE DESCRIPTION 

FREQ 

GEN 

MASS 

DEFL 

FREQ 

GEN 

MASS 

DEFL 

WING ENGINE YAW 

46 

77 

- 

55 

- 

- 

FIRST WING BENDING 

47 

78 

136-143 

56 

86 

166-173 

AFT ENGINE YAW 

48 

79 

- 

- 

87 

- 

HORIZONTALSTABILIZER BENDING 

49 

80 

144-149 

57 

88 

174-180 

VERTICAL STABILIZER BENDING 

50 

81 

150-157 

59 

89 

189-196 

AFT FUSELAGE TORSION/ 
WING ENGINE PITCH 

51 

82 

158-165 

58 

- 

181-188 

AFT FUSELAGE LATERAL BENDING 

52 

83 

- 

- 

- 

- 

SECOND WING BENDING 

53 

84 

■* 

60 

90 

- 

SECOND VERTICAL STABILIZER 
BENDING 

i ■ 

54 

85 

_ 



- 

- 

- 
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FIGURE 25. MODE SHAPES AND MODE LINES 



FIGURE 26. MODE SHAPES AND MODE LINES 
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FIGURE 33. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 32. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 34. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 35. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 37. EVOLUTION OF FREQUENCY RATIO 



FIGURE 36. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 38. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 41. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 40. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 42. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 45. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 46. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 47. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 49. EVOLUTION OF FREQUENCY RATIO 



FIGURE 48. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 50. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 51. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 53. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 52. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 54. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 55. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 57. EVOLUTION OF FREQUENCY RATIO 


. r .,. : 
4- i- 


■!---!• .L.l.. 


nding ■ }. j: . t . ;.. 

■c I'uuW ' •/!.: ! 4 I • L .4 


M : I ; J 

Itt.WlNlG.GEI 


AHTISI^tTNIC | FULL FUEL i ; | J | 

I ! .: m;..! j. ; :/ ±2 . 

<5 : <j>: | : 


4 ' I ' : 1 M 4 1- I i rl: 


4:ii;t[no4 


J .j WING STIFFNESS ■ 


I. I 

|-r' + 


: r; 


4 - 4 1 4> 


L 4. L 1.4 

\±L 


i 0,2 

Ltr.ru 


WING PILON STIFFNESS : \ • ... ... ... I. J .j I j j: 

ENGINE WEIGHT AND FUSELAGE STIFFNESS ! j H L ..{ f- iL i . L 

w ng stiffness and aft ptlon Stiffness . ; I... 1 

W NG ENGINE WEIGHT AND WING PYLON STIFFNESS ,|_J j _4_ „ ..Li 

I W NG STIFFNESS I \ ! !' 44. 

... I WING STIFFNESS AND FUSELAGE STIFFNESS .j ! j- . ..lilii 

I WING STIFFNESS I \ . j , 1 4 

| ... . _ _ WING STIFFNESS, i | _ . ; 

I WING PILON STIFFNESS i I 1 ! I , 

I _ 4 .. ...I... [ AFT ENGINE WEIGHT AND PILON STIFFNES; 

1 li lit; ■ j :l;i vjrtft II 

• 4 ' 6.8 ]0 12 : If* 'l<j 16 ; 2f> ' ^2 ! 24 ; 26 • 2p ‘ ^0 j 3? i 34 ^ 44 

. . .. [Q ± 


FIGURE 56. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 58. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 59. EVOLUTION OF FREQUENCY RATIO 
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FIGURE 6T. EVOLUTION OF GENERALIZED IWASS RATIO 
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figure 62. EVOLUTION OF GENERALIZED MASS RATIO 
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FIGURE 63. EVOLUTION OF GENERALIZED MASS RATIO 


FIGURE 64. EVOLUTION OF GENERALIZED MASS RATIO 




FIGURE 65. EVOLUTION OF GENERALIZED MASS RATIO 


FIGURE 66. EVOLUTION OF GENERALIZED MASS RATIO 




FIGURE 69. EVOLUTION OF GENERALIZED MASS RATIO 



FIGURE 70. EVOLUTION OF GENERALIZED MASS RATIO 
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FIGURE 73. EVOLUTION OF GENERALIZED MASS RATIO 


n * t i ; ! 

r '- 1 


STfctflk! : L (VIS. I 


t rrpr 


, r .... i r - 

M'i .1-4 — l 1 — !. - 

1M ; i i n . 


l<W ‘ , 


■I i i -t M^ T 

li -1 -h-l'ni. 


m ilir Hn rr m rr 4 fin|iMg 

, if Mi RAl I II D .j, .1 J. , , 1 . . _ j. 

m ass i | ' i r i . ' L- , . • • 

I MTIO I . 'll , , , 

<«v«y r •* I 1 • ■ ;- < ! t '| T ■ " t 

i , 1 .» ! . -..I .j, 

I NIG PYLON STITFHtSS | 

■ 4 ; l I ; I CKUC WtlCtfT^HD fuSELAGt STIFRESS 1 \..J r 

I | | l°'°. V HG STIFFNESS ttT P^ON $T I Ft JESS ' | 

. U NG ENGINE VtlGHTMD WIMlPrLOi STIIFMEIS ; . 

; . ; v *& s rimecss ■' / , 

L I I . _i WING surf N<SS AND ruSELAdE ST ^FFXEtS[ 

. ; i . wing srirmiss I 

] _! I ■ 2 _i “WS STIFFNESS J 

f ■ I ” 1 1 1 | MING PTLDN sjlFFKESS 

1 ! I 1 1 •[.I AFT ’ EKQ] Mt WEIGHT W 

! i l M l cmiir<c»i 


mt i « 


"L FVllirv^ J'irrnuj . . | . . *. - r _ i. - . , 

I *; viFFmssffetJ ttr p&qk !HF(i£Ss I j I 
y pic CXOIIIt WtieKTABD WIHCIPriOH STIIFHEIS , |. 

( v « inir*Fi5 I | I 1 I , 

H1K SIlFFIlCSS «D FUSll /(4 5l(FF>Eksj 

“Ip‘ I 


ymp Pruw siiFFutsj 1 i I . ; 

. I AFT.EKdht utiOHT wo tarn slir/yfsl 


.1 I • 

M ! , I -i 


| 1 ' CEpUF|CATl|N (Gp| 

0 l| . lift lil !<> 2? ' J* 21 jp 1 :o 3 I If | i I 

WH'jrrl ill tii-l J 1. Ha 11 h I ! 


FIGURE 72. EVOLUTION OF GENERALIZED MASS RATIO 
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FIGURE 74. EVOLUTION OF GENERALIZED MASS RATIO 



FIGURE 77. EVOLUTION OF GENERALIZED MASS RATIO 
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FIGURE 79. EVOLUTION OF GENERALIZED MASS RATIO 



FIGURE 81. EVOLUTION OF GENERALIZED MASS RATIO 



FIGURE 80. EVOLUTION OF GENERALIZED MASS RATIO 



FIGURE 82. EVOLUTION OF GENERALIZED MASS RATIO 



FIGURE 83. EVOLUTION OF GENERALIZED MASS RATIO 
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FIGURE 85. EVOLUTION OF GENERALIZED MASS RATIO 
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FIGURE 84. EVOLUTION OF GENERALIZED MASS RATIO 
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FIGURE 86. EVOLUTION OF GENERALIZED MASS RATIO 
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FIGURE 105. EVOLUTION OF DISPLACEMENT RATIO h 
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FIGURE 109. EVOLUTION OF DISPLACEMENT RATIO a 
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FIGURE 122. EVOLUTION OF DISPLACEMENT RATIO h 
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FIGURE 125 . EVOLUTION OF DISPLACEMENT RATIO h 
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FIGURE 133. EVOLUTION OF DISPLACEMENT RATIO h 



FIGURE 132. EVOLUTION OF DISPLACEMENT RATIO a 
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FIGURE 134. EVOLUTION OF DISPLACEMENT RATIO h 
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FIGURE 135. EVOLUTION OF DISPLACEMENT RATIO a 
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FIGURE 137. EVOLUTION OF DISPLACEMENT RATIO 4> 
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FIGURE 136. EVOLUTION OF DISPLACEMENT RATIO 0 



FIGURE 138. EVOLUTION OF DISPLACEMENT RATIO h 




FIGURE 139. EVOLUTION OF DISPLACEMENT RATIO a 



FIGURE 141. EVOLUTION OF DISPLACEMENT RATIO a 




FIGURE 142. EVOLUTION OF DISPLACEMENT RATIOS 
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FIGURE 151. EVOLUTION OF DISPLACEMENT RATIO <// 



FIGURE 153. EVOLUTION OF DISPLACEMENT RATIO a 
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FIGURE 154. EVOLUTION OF DISPLACEMENT RATIO h 




FIGURE 155. EVOLUTION OF DISPLACEMENT RATIO a 



FIGURE 157. EVOLUTION OF DISPLACEMENT RATIO * 
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FIGURE 156. EVOLUTION OF DISPLACEMENT RATIO 9 . 



FIGURE 158. EVOLUTION OF DISPLACEMENT RATIO 0 
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FIGURE 161. EVOLUTION OF DISPLACEMENT RATIO a 
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FIGURE 160. EVOLUTION OF DISPLACEMENT RATIO h 
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FIGURE 162. EVOLUTION OF DISPLACEMENT RATIO h 
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FIGURE 176. EVOLUTION OF DISPLACEMENT RATIO h 


!■ ! • • P i • -i 

itoRIll. <T« flip Wf^TtrAW DLfLECTipN (W) 

j ; | : 

• i I : ; : 

HDRIZ star, boding :3.oo|h*.|_- j 

. I 

i. 7[- |j 

8NTISWETRIC [ jFULL |FVtL; j 


i- i 1. r i 


•1.00 • <& (|)©^®- j- --0-1-$— ^ 


DISPLACEMENT ...i .B 
RATIO | 

(h/h.) 1 | 


!ESS i . J..,..- ] _ ,_ij jjl J.. 


I ’ ! • : m 

■ , 4- - !•••••• 4 


i J . 6 WIHG STIFFHE5S : I I . [ \ ! l^„l I I : | ! 

j ! " j ! I WING PYLON STIFFNESS i : i 1 i !“ [7 T " i j T" 

j L ..i .... ENGINE WEIGHT AND FUSELAGE STIFFNESS ! .. j . . .Li. ..| 

( I ! . ! WING STIFFNESS AND AFT PYLON STIFFNESS : . . .! i j' " 

i . .L. ■ ..1. ..! . ! .4 .. U NG ENGINE WEIGHT AND WING. PYLON STIFFNESS .. S : j j . L. : . 

j . . | . w ng stiffness > ■ i i* ! l ' : .; 

| -_J _ J ' WING STIFFNESS AND FUSELAa STIFFNESS . ’ .1. J-L. 

I ... I ! ! MING STIFFNESS ! . . . I , : ] 1 - ! -i . I . . J-I 

-[-■ lJ ...„ L 2 L I WING STIFFNESS I ^ j 1 ......l: ~ . ... i » 

' T -1 ! * 7 : WING RYL^N STIFFNESS ! ' I*. 

L ; 1 . J. ■ ! . . . : AFT ENGINE WEIGHT ANO PYLON STIFFNfS! 

" ! 'll CERTIFICATION (GVT) i ■ 

■ .j- -i - : i-i -i-i - 1 I . ! ; m i 

| • ■■■"- 0 - j |j f f 8 10 II 1|4 16 18 20 22 24 26 28 • 3!: ; M If yn 

: , ! I.. ■ |Ti ! , ! l — l i— .Tw'A.MNTHS.-,Mj 1 [' — i— J : I I ■ j ..I.... 


. l h . p . : . .l i ....!.: . . • .~tzj .. i... . 1 n mn 

FIGURE 178. EVOLUTION OF DISPLACEMENT RATIO h 



Pin rrp ; r : : • I; i.'r;:u |J..i w;x:r :r i:; r i- ! j 7 ; i ; 

" 1 j j • I- ; f r M :■ i iopi* ' ! : I i : i ' 1 ■ ; - H i- 

iTjf 


-rrf7 i' 

r T IT 7 


-- ■ - I : ■ j-' 

f4 oisplacenent I-4 ■* i ~i - , 

> AATIO j. I , I . ■ 

r: <'V I- - ; 


j ':«n»MnK! i IIU ML ’ : • . ; J j ; ■ • 

r: ' T i-T-.! ;-i- ; i - 

r -- - : ---i-H- f J r r 

I ; • • • j ] : 1 ■ ■ i 

: r jZ~:p;~ip' '""“I'ijr ft" j" 



| _ .j HIM STIFTHC 55 !_ j } j }; | r T" -| j |~ ! J 

T HIM pvLm~ STIFFNESS I !. | l . \ | . "ITT 

j. . . J... ' EMC INC HEIGHT iAM fUSEUSE STIFFNESS L .;. '■ I •_ 

j ! HIM STIFFNESS WO AFT PYLON STIFFNESS { 

-L. .4 HIM ENGINE HEIGHT ANO HIM PTLOR STIFFNESS 

j . HIM STIFFNESS = I ! I A, 
4-. J . < ><L. L. . HIM STIFFNESS. AM) FU5ELA8LSufFJ 


STIFFNESS. MO FUSELAflE.S ufFNESS . L_ 

T«Sfe.:l! T i 


xuiiu r. 


' HIM FYL'ON STIFF'S! I' ll: !.. 

» 1. _ ... . I - ; AFT'EMIM WEIGHT AM #lfL0N! STIFTNPS! 

I I i CERTIFICATION (GET) I 

- 44 ■ ; — I T I 4- ! — I — I — J-4M-4 

-- 4 — 4 - Itw^Ncimai 7 r ~~t 1 -J ; ]■ 7 17- 


FIGURE 179. EVOLUTION OF DISPLACEMENT RATIO £ 

■ «•». X-, ;. t r ft-rj-j ! • • 


'• j 5 j : i : < 

: ! , ! J 11.2 


! i 


D1SPLAC£I€HT 
RATIO 
( / J 


11.2 j 

: I 

I ■: • 

j,:. ... 

.i-i. . 

i .. i: 


, I ! ! i l I ! i ■ II 

ri :• Iri” -hr- ; ' ! * : 

j . . : ! nfawt hose nou amle [*) . 

' AFT FUS :UU£; TTWIOH/ltlNG fuel* PI^CH. 3.36 lit I 

’ ! (wuwermt • j Mu. ml: j ■ j : 

1 |."T”®j'r i: 7 j"7'l j f ! I ; ! 

i "i i t!ttt : r i- : ;■ 


., WING STIFFNESS ! jJ.J7Jj.4_.: ! .. ]._. LJJ..:. 1 

mk mo* STIFFNESS : T : [ I . | I 

ENGINE HEIGHT AND FUSELAGE STIFFNESS ! . j ...i L... 

KING STIFFNESS MD AFT FTLON STIFFNESS : I . [ 

.4 WING ENGINE HEIGHT AND UING FTLON STIFFNESS . 

| WING STIFFNESS ! 1 

! WING STIFFNESS ANO FUSELAGE STIFFNESS . |. 


; | j j . I ..4 


4;i 



HIM PYLON STIFFNESS 1 ! ' j 

AFT ENGINE HEIGHT AM PYLON ST I FFMS! 
•It CERTIFICATION (GVT) i 


I..L.L.U ,j.j i.:i4 


^ 10 12 14 If 20 22 ?4 24 2S 30 32 

7— ! L !tine!% months : _j — j — 1 — i_«J I 


FIGURE 181. EVOLUTION OF DISPLACEMENT RATIO d 


vO 
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SECTION 4 


SOME POSSIBLE APPLICATIONS 

As noted previously, this report is intended to provide a data base for the 
control system designer so that he can estimate the likelihood that an 
advanced system will provide its intended function at first flight in the 
presence of uncertainties in plant and control system dynamics. The specific 
application of these data will depend on the designer's favorite analytical 
tools and the general design approach used by a manufacturer. This section 
will suggest a few ways of using these data. As will be apparent, the assess- 
ment of stability margins will be reasonably straight-forward, but the 
estimation of flying qualities is not quite as easy. Reference 4 gives some 
results of recent transport simulator handling qualities tests. It is 
apparent that conflict exists between various popular handling qualities 
measures, so the system designer will need to use several yardsticks to 
measure the effects of these reported parameter variations. 

CLASSICAL SERVO ANALYSIS APPLICATIONS 

Normally the system designer uses classical analysis methods to a first-cut 
definition of a system based on nominal aerodynamics and an exact specifica- 
tion of control dynamics. A second iteration then considers the change in 
closed-loop dynamics and system performance due to plant and controller 
parameter variation. The classical analyst depends on the visibility of 
system characteristics via root-locus. Bode, and time history plots to aid 
him in the system design, so it would be desirable to reflect the effects of 
parameter variability graphically. 

If we consider a simple pitch axis problem, the aircraft nominal open-loop 
dynamics might appear as shown in Figure 197. 

The tuck/subsidence, short-period, actuator and first structural modes 
(nominal) are shown along with the pitch to elevator numerator. The control 
system designer may be interested in an attitude-hold system for this air- 
craft and would then make the usual rate plus position closure as shown in 
Figure 198. 
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FIGURE 197. TYPICAL LONGITUDINAL ROOT LOCATION 


The nominal closed-loop system is well behaved and performs to specification, 
but the effects of variations in the system parameters are not obvious. 


Past work has established the mathematics for assessing the effects of 
uncertainties on closed-loop characteristics (e„g., Reference 8) but these 
methods have been so tedious for higher order systems that they have not been 
generally applied. Computer capacity and speed are now sufficient to allow 
sensitivity methods such as those developed in reference 8 to be programmed. 

A Monte Carlo parameter variational approach could be used to establish the 
"envelopes" of open-loop and closed-loop system parameters for expected plant 
parameter variations as sketched in Figure 199. 
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FIGURE 198. TYPICAL LONGITUDINAL 
SYSTEM ROOT LOCUS 



FIGURE 199. CLOSED-LOOP ROOT 
VARIATION 


s 


The designer's task will be to design the system such that e mi - n < ? < 
etc„, to satisfy the basic flying qualities and structural criteria. Where 
the modes being stabilized are potentially "flight critical" (e.g., the first- 
order divergence or second-order low-damped modes shown in Figure 199 could 
possibly result in unacceptable characteristics), the designer will take all 
possible assurances that the total pilot/vehicle system will perform 
reasonably well before the first flight opportunity arrives. 

STATE-SPACE SYNTHESIS METHODS 

In an initial design exercise, state-space methods are now often used in the 
form of linear-quadratic-Gaussian synthesis procedures for quickly identifying 
candidate control systems (References 5 and 6). With estimates of plant 
uncertainties available, state-space methods can be employed to determine 
eigenvalue sensitivities in addition to time response and frequency response 
sensitivities already discussed. Sensitivities may be determined for 
individual as well as statistically selected combinations of parameter 
variations. Alternative control system configurations which yield approxi- 
mately equivalent dynamic performance characteristics are easily synthesized 
using state-space methods. The sensitivity of these systems to plant param- 
eter variations offer criteria for selection of a particular configuration. 

Eigenvalue sensitivities can be obtained from a pole placement synthesis 
program which employs a first order perturbation approximation to relate 
changes in the closed-loop system coefficient matrix. A, to changes in the 
system eigenvalues , A. 

The matrix equation employed is: 

D Aa = AA 

Where D is the eigenvalue sensitivity matrix, Aa is a vector consisting of the 
system parameter variations, and a a is a vector consisting of the eigenvalue 
changes. The sensitivities of the eigenvalues to parameter variations from 
the nominal values may be determined merely by examining the matrix, D. D 
would be calculated for various control configurations, and would also be 
recalculated for cases where the parameters are varied from the nominal 
values. '(Note '.that .a a consists of the non-zero elements of the matrix aA). 
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The off-nominal sensitivity matrix is obtained by calculating D using the 
matrix (A + aA) rather than the nominal matrix A.) 

PARAMETER IDENTIFICATION AND ADAPTIVE CONTROL APPROACH 

Yet another system design philosophy could be to conduct an exercise that 
shows the ai rcraf t/control system has closed-loop dynamics meeting the 
designer's specification for the entire range of plant dynamics described by 
the plant parameters with their expected variation. Such a system might be 
required to converge within a fixed period of time (which would be much less 
than the time for control problems to develop) from any and all corners of the 
open-loop dynamic envelope of Figure 199, and at all flight conditions. The 
approach would represent a significant departure from the previous two 
approaches where control system gains are assumed fixed at a specific flight 
condition. 

Such a system design would employ parameter identification techniques to 
determine plant parameters and to update the control laws in response to on- 
line measurements of the plant states. A comprehensive bibliography on the 
system identification problem is contained in Reference 7. 

POSSIBLE STATISTICAL USAGE OF STRUCTURAL DYNAMIC RESULTS 

The previously described study utilizes only one airplane as a basis for 
parameter variability determination. The lack of data for a large number of 
airplanes is compensated to a great extent by the consideration of parameter 
variability within a reasonably large group of different flexible modes. 

For each modal parameter, a measure of variability can be established by 
bringing the results together from each of the modes. Let r-jj(t) denote the 
value at time t of any of the parameter ratios expressed in the plots. The 
subscripts indicate that the result for the i-th parameter is being considered 
for the j-th flexible mode and r^(t) -*■ 1 as t -*■ T, where T denotes the 
durations to certification. If b-j(t) signifies the average of the logarithms 
of the i-th parameter values over the modes (i.e., over j), then 

b i (t) ■ H | 
j * . 
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M is the number of modes presented in the plots in a defined category (e.g. 
symmetric, zero fuel) and b. (t) is the bias tendency of the log of the i-th 
parameter at time t in the design cycle. Then the variability can be measured 
by v^(t), where 

V 2 (t) = (jjprj * Clog r. ^(t) - b.(t)] 2 

Having the bias and variability estimates for each parameter, combinations of 
parameter magnitudes can be selected for investigating the control system 
design. In studying the airplane performance to control and external inputs, 
the parameters used for the airplane can be selected as follows: If p . . (t) 

* J 

is the value determined at time t for the i-th parameter of the j-th mode, 
then a set of nominal parameter values n..(t) can be obtained as 

* J 

★ 

log n.jft) = log P-jj(t) - b..(t); i = 1, .... N 

★ 

j = 1 , .... M 

A set of extreme parameter values e..(t) are given by 

■ J 

log e ij .(t) = log p ij .(t) - b.(t) + 3 v.(t); i = 1, ..., N* 

j = 1 M* 

(The number of airplane flexible modes M* used for the control system design 
may be different than the number M defined above. N* is the number of 
parameters pertinent to the M* modes). 

The determination of parameter combinations from the P = M* • N* parameters 
can be done as follows: In one combination all parameters can be set at their 

nominal values. For the other 2P combinations, (P - 1 ) of the parameters 
are set at their nominal values while the remaining parameter is set at each 
of its extremes. This simple scheme gives 2P + 1 parameter combinations. 

In a previous discussion, it is indicated that a Monte Carlo procedure can be 
used in selecting the parameters. In this approach it can be assumed that the 
log r^Ct) have independent normal distributions with mean b.(t) and standard 
deviation v,.(t). From a table of random numbers, values of log r.j(t) are 
drawn for each of the P parameters. Now p. .(t) is the value of the parameter 

* J 
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determined at time t in the design cycle. Let p^(t) be a randomly deter- 
mined value for the sape parameter. This ..value can be found as follows from the 
randomly selected r. - (t) value: 

Pij (t5 = n u (t) • r ij (t) 

where n.^^t) is previously defined. The number of combinations required for 
a Monte Carlo approach will at least be as large as that used-, in the above 
simple scheme. 


The usages of the plotted results that have been discussed above can be 
further clarified by considering the following airplane openloop differential 
equations for the flexible modes: The equation for. the j-th mode, in symbolic 
form, is - 




(t) 


+ (4Tr 2 6 aj ffj > j + q a (t)] = . T Fg(t) 0 Bj . 


(s • = 1 for a = j; = 0 otherwise) j = 1 , ... , M* 

J 


This second order equation can be converted into two first order equations for 
the state-space formulation. This transformation is straight-forward and, will 
not be done here. The parameters expressed in the above equations represent 
the following for the j-th mode: 


m. - generalized mass 

J 

t. - critical damping factor 

J 

f . - natural frequency 

J 


0 


0 ^ - modal coordinate at location 


no 



C j«* C j«* C j« - generalized unsteady aerodynamic coefficient matrix 
elements associated, respectively, with aerodynamic 
j , forces produced by accelerations, velocities and 
. displacements. 

The time-dependent quantity q^t) is the generalized displacement for thea-th 
mode. q o {t) and q tt (t) are the associated first and second time derivatives. 
Fg(t) Is an external time-varying force acting on the mode shapes at the 
number 6 modal coordinate. 

The number. A, of modes over which the left-hand side summation Is conducted 
equals the number of flexible plus rigid body inodes being considered. The 
number of modal coordinates B for the right-hand side summation Is given by the 
size of the modal eigenvectors. The number of parameters N* necessary to 
define each modal equation Is at most 3 + 3A + B. (A modal coordinate 0 6 . 
is not needed if no force is associated with it). 


The N* parameters for the j-th mode can be placed in the j-th column of a 
matrix array. The 1-th row element of this column is the parameter quantity 
p^i previously defined. Since the parameters pertain to a particular time t 
in the design, they are designated as P^j(t). No data has been presented In 
this report for the parameter P^j(t) = Cj to calculate Its bias or variability. 
Its value selected for any control system study will then remain fixed at its 
determined value. A similar statement applies to each of the coefficients 
C ja* Cj, and C^. Thus of the possible set of N* parameters for the j-th 
modal equation, the results In the given plots will provide bias and varfabllty 
information for at most 2 + B of the parameters. 


CONCLUSIONS AND RECOMMENDATIONS 


Based upon the uncertainty investigation over the design period of a single 
large transport aircraft the following conclusions are made: 


1. The largest variation In modal parameters occurs during the early design 
phase, reflecting major design changes. 

2. Analyses performed after the early design phase show less variation in 
modal parameters, reflecting minor modification to the airplane design data. 


Ill 


3. No dlscernable trend was seen In the parameter uncertainties between fuel 
condition or airplane symmetry. 

4. The modal generalized mass and displacement parameters show much larger 
variations than do the modal frequency parameters and are of limited use 
except as input to further response analyses. 

5. The magnitude of the input modal displacement ratio shows little or no 
correlation to the magnitude of the output response acceleration ratio. 

6. By the time the basic design characteristics have stabilized, the struc- 
tural dynamic modal characteristics can be predicted adequately for use 
in the control system design. 

7. The generalized mass and modal amplitudes, by themselves, do not provide 
insight to the airplane flight response characteristics. 

The modal parameter uncertainty results presented herein represent a starting 

point for further studies. Recommendations are as follows: 

1. The preferred location of sensors is often at nodal or antinodal points. 
The uncertainty of these points should be determined to enable optimum 
sensor placement at an early part of the design cycle. 

2. The effect of variations in airplane (i.e. plant) parameters on airplane 
flight control should be evaluated in a follow-on phase. The parameters 
whose variations produce the most pronounced effects should be identified. 
Control system designs should then be studied to determine the extent to 
which these sensitivities can be reduced. The proper positioning of 
sensors and control surfaces would be part of this study along with the 
design of the control system. 

3. The effect of unsteady aerodynamics uncertainties should be investigated 
for the purpose of not only improving control system functions but for 
the immediate problem of flutter prevention. 
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APPENDIX A 


-%x 

BODY AXES EQUATIONS OF MOTION 

Longitudinal Set: 

u = _ W Q q - g9 cos 0 Q + X u u + X q q + fy* + X*w + X 6{; 6 e 

* = U 0 q - g© sin 9 q + Z u u + Z q q + Z w w + + Z 6e 6 e 

3 = M y u + M q q + M w w + + M 6e s e 

Lateral -Directional Set: 
v 

P 
r 


= - U Q r + W Q p + g ip sin 0 Q + G 4> cos 9 Q + Y r r + Y y v + Y p p 


+ Yx 6 + Yx <5 + Yx 6 

6 a a 6 r r 6 sp sp 


= r + Lr + Lv + Lp + Lx6 + Lx 5 + Lx 6 

I r v p K 6 a a 6 r r °sp sp 


Pr + V + V + V + N 6 a 6 a + N «rV + N «sp 6 


sp sp 
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APPENDIX B 


DEFINITION OF COEFFICIENTS OF EQUATIONS OF MOTION 
(DIMENSIONAL STABILITY DERIVATIVES) 


W * a_U 
0 0 0 


pUSb r 
4m C y„ 


X = — — f-C - c 
A u m 1 L D u c D 

y - pUScT r 
A q ' 4m 


\ -'*£• Cy 


p 4m, V 


X w = 2m S (C L ’ C dJ 


Y. 


2m V, 


y _ P SC p 

A w “ " 4m 


t . -ayfic 
V 2" h 


V _ PITS . 
\ ' 2" \ 


I - PUSb 1 
r " 41.. 


4I x C * 


1 = £y§. (r _ r ) 

u m { l L u V 

7 1 pUSc r 

q 4m u L q 

7 = p U ? (_q _ r ) 

4 w 2m 1 l L„ V 


L c 

2I X % 

I = .&ySb? c 

" «* S 

, _ pU 2 sb r 

\ ' 2I X V 


7 — P p 

w " " 4m 


aU 2 Sb r 

2I X ‘5. 


Z .J^c 
°e 2m C L, 


N * pusb r 

N v " ^17 C n 


M = ^ (C + C) 
u I y m u m 


pUSb 2 r 
C n r 
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M q 

"w 


pUSc 2 r 
' 4r y S 
_ pUSc r 


2I y m 


"w 


pSc 2 r 

4Iy X 


pU 2 Sc r 

2I y H 


PUSb £ r 

N P -^7 C n p 

n, >-^ a -c„ 

5 r 2I z 


pU 2 Sb „ 

Nj - — of 6_ 

6 a 2I z n 6. 


N = pU ^ C 

V n 21 V 

S P z 6 sp 

. _ pU 2 Sb r 

L S 21 a. 
sp x 6 


sp 


V - r 

p “ 2m V 

SP 6 sp 
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